clear all
set printcolor gs1, permanently
set more off
global dirout ""
global dirgeolocalized ""
global dircenso ""
global ENIF = ""


use "$dirgeolocalized/users_geolocalized.dta", clear
gen users=1

collapse (sum) users count, by(ID)
rename ID_origin ID

merge 1:1 ID using "$dircenso/censo_CDMX_geo.dta", nogenerate keep(match) keepusing(POBTOT GRAPROES VPH_AUTOM VPH_CEL VPH_INTER VIVTOT )


foreach var in VPH_AUTOM VPH_CEL VPH_INTER {
gen s`var'=`var'/VIVTOT
}

keep if VIVTOT>0
keep if POBTOT>0


gen share_users=users/POBTOT 
replace share_users=1 if share_users>1

		binscatter share_users GRAPROES if GRAPROES>5 & GRAPROES<16, ///
		nq(100) linetype(none) lcolor(navy) mcolor(ebblue) xtitle(Average years of schooling) ///
		ytitle(Share of users)savegraph("$dirout/FigureC2a.eps") replace  
		
		binscatter share_users sVPH_INTER [w=VIVTOT], ///
		nq(100) linetype(none) lcolor(navy) mcolor(ebblue) ///
		xtitle("Share of homes with internet") ///
		ytitle(Share of users) savegraph("$dirout/FigureC2b.eps") replace
		
		binscatter share_users sVPH_CEL [w=VIVTOT], ///
		nq(100) linetype(none) lcolor(navy) mcolor(ebblue) ///
		xtitle("Share of homes with cell phone") ///
		ytitle(Share of users) savegraph("$dirout/FigureC2c.eps") replace
		
		binscatter share_users sVPH_AUTOM [w=VIVTOT], ///
		nq(100) linetype(none) lcolor(navy) mcolor(ebblue) ///
		xtitle("Share of homes with car") ///
		ytitle(Share of users) savegraph("$dirout/FigureC2d.eps") replace
		
 	

cd "$ENIF"

use  "TVivienda.dta", clear
keep UPM VIV_SEL ENT MUN 
save "TVivienda_2015.dta", replace

use  "tmodulo1.dta", clear
keep UPM VIV_SEL HOGAR SEXO R_SEL EDAD NIV P3_1 P3_2 P3_5 P4_8 P3_8_1 P3_8_2 P5_21 P3_11 P3_12 FACTOR 
joinby UPM VIV_SEL using "TVivienda_2015.dta", _merge(matched)   
tab matched
drop matched
gen clave_municipio=ENT+MUN
rename ENT clave_estado
drop MUN UPM VIV_SEL HOGAR R_SEL
destring SEXO P3_1 P3_2 NIV EDAD P3_5 P3_8_1 P3_8_2 P3_11 P3_12 P4_8 P5_21 FACTOR clave_municipio clave_estado, replace 

rename SEXO sexo
replace sexo=0 if sexo==1
replace sexo=1 if sexo==2
label define sexol 0 "Men" 1 "Women" 
label values sexo sexol

rename P3_1 relacion
label variable relacion "¿Qué es usted del jefe del hogar? jefe=1, esposo=2, hijo=3, nieto=4, otro=5, sin=6"
 
rename P3_5 trabajo
label variable trabajo "trabajo=1, sintrabajo=2, buscotrabajo=3, estudiante=4, hogar=5, jubilado=6, incapacitado=7, no=8"
gen trabajo_bins=.
replace trabajo_bins=0 if trabajo==5
replace trabajo_bins=1 if trabajo==6
replace trabajo_bins=2 if trabajo==4
replace trabajo_bins=3 if trabajo==2 | trabajo==3 | trabajo==7 | trabajo==8 
replace trabajo_bins=4 if trabajo==1
label define trabajo_binsl 0 "Work at Home" 1 "Retired" 2 "Student" 3 "Unemployed" 4 "Employed" 
label values trabajo_bins trabajo_binsl

rename P3_2 estado_civil
label variable estado_civil "union libre=1, separado=2, divorciado=3, viudo=4, casado=5, soltero=6"
gen casado=0
replace casado=1 if estado_civil==5
label define casadol 0 "Single" 1 "Married" 
label values casado casadol

rename NIV educacion
label variable educacion "no=0, kinder=1, prim=2, secun=3, tecysecun=4, normal=5, prepa=6, tecyprepa=7, lic=8, masdelic=9"
gen educacion_bins=0
replace educacion_bins=1 if educacion==2
replace educacion_bins=2 if educacion==3 | educacion==4 
replace educacion_bins=3 if educacion==5 | educacion==6 | educacion==7
replace educacion_bins=4 if educacion==8 
replace educacion_bins=5 if educacion==9
label define educacion_binsl 0 "None" 1 "Primary" 2 "Secondary" 3 "HS" 4 "BA" 5 "BA+"
label values educacion_bins educacion_binsl

rename EDAD edad
replace edad=. if edad==99
label variable edad "18-70"
gen edad_bin=1 if edad<=34
replace edad_bin=2 if edad>34 & edad<=44
replace edad_bin=3 if edad>44 & edad<=54
replace edad_bin=4 if edad>54 & edad<=64
replace edad_bin=5 if edad>=65
label define edad_binl 1 "<35" 2 "35-44" 3 "45-54" 4 "55-64" 5 "65>" 
label values edad_bin edad_binl

rename P3_8_1 ingreso
label variable ingreso "¿Cuánto gana o recibe usted por trabajar (su actividad? 999999=NA)"
replace ingreso=. if ingreso==99999

rename P3_8_2 ingreso_freq
label variable ingreso_freq "¿Cada cuándo le pagan? 1=semana, 2=quincena, 3=mes, 4=año, b=blanco"

rename P3_11 celular
label variable celular "¿Usted tiene celular? 1=si, 2=no"
replace celular=0 if celular==2

replace ingreso=ingreso*4 if ingreso_freq==1
replace ingreso=ingreso*2 if ingreso_freq==2
replace ingreso=ingreso/12 if ingreso_freq==4
drop ingreso_freq

rename P3_12 celular_aplicaciones
label variable celular_aplicaciones "¿Su celular permite descargar aplicaciones como whatsapp, facebook o twitter? 1=Si, 2=No, 3=No sabe, b=blanco"
replace celular_aplicaciones=. if celular_aplicaciones==9
replace celular_aplicaciones=0 if celular_aplicaciones==2

rename P5_21 utilizo_debito
label variable utilizo_debito " En los últimos tres meses ¿usted utilizó su(s) tarjeta(s) de débito? 1=Si, 2=No"
replace utilizo_debito=0 if utilizo_debito==2
gen households=1
gen smart_and_debit=celular_aplicaciones*utilizo_debito

save temp.dta, replace

xtile pct=ingreso [aw=FACTOR], nq(10)
collapse (mean) celular celular_aplicaciones utilizo_debito smart_and_debit [pw=FACTOR], by(pct)
drop if pct==.


gen utilizo_debito2=utilizo_debito-smart_and_debit
gen celular2=celular-celular_aplicaciones
replace utilizo_debito2=0 if utilizo_debito2<0


gen celular_aplicaciones2=celular_aplicaciones-smart_and_debit
graph bar smart_and_debit celular_aplicaciones2, over(pct) stack ///
bar(1, color(cranberry) lcolor(black)) bar(2, color(midblue) lcolor(black))  ///
legend(label(1 "Debit & Smartphone" ) label(2 "Smartphone")) ///
ytitle("Share of households", height(5)) ///
graphregion(color(white)) bgcolor(white) 
graph export "$dirout/FigureE1d.eps", replace 

use temp.dta, clear
collapse (mean) celular celular_aplicaciones utilizo_debito smart_and_debit [pw=FACTOR], by(edad_bin)

gen celular_aplicaciones2=celular_aplicaciones-smart_and_debit
graph bar smart_and_debit celular_aplicaciones2, over(edad_bin) stack ///
bar(1, color(cranberry) lcolor(black)) bar(2, color(midblue) lcolor(black))  ///
legend(label(1 "Debit & Smartphone" ) label(2 "Smartphone")) ///
ytitle("Share of households", height(5)) ///
graphregion(color(white)) bgcolor(white) 
graph export "$dirout/FigureE1f.eps", replace 


use temp.dta, clear
collapse (mean) celular celular_aplicaciones utilizo_debito smart_and_debit [pw=FACTOR], by(educacion_bins)

gen celular_aplicaciones2=celular_aplicaciones-smart_and_debit
graph bar smart_and_debit celular_aplicaciones2, over(educacion_bins) stack ///
bar(1, color(cranberry) lcolor(black)) bar(2, color(midblue) lcolor(black))  ///
legend(label(1 "Debit & Smartphone" ) label(2 "Smartphone")) ///
ytitle("Share of households", height(5)) ///
graphregion(color(white)) bgcolor(white) 
graph export "$dirout/FigureE1e.eps", replace 


use temp.dta, clear
collapse (mean) celular celular_aplicaciones utilizo_debito smart_and_debit [pw=FACTOR], by(sexo)

gen celular_aplicaciones2=celular_aplicaciones-smart_and_debit
graph bar smart_and_debit celular_aplicaciones2, over(sexo) stack ///
bar(1, color(cranberry) lcolor(black)) bar(2, color(midblue) lcolor(black))  ///
legend(label(1 "Debit & Smartphone" ) label(2 "Smartphone")) ///
ytitle("Share of households", height(5)) ///
graphregion(color(white)) bgcolor(white) 
graph export "$dirout/FigureE1b.eps", replace 


use temp.dta, clear
collapse (mean) celular celular_aplicaciones utilizo_debito smart_and_debit [pw=FACTOR], by(casado)

gen celular_aplicaciones2=celular_aplicaciones-smart_and_debit
graph bar smart_and_debit celular_aplicaciones2, over(casado) stack ///
bar(1, color(cranberry) lcolor(black)) bar(2, color(midblue) lcolor(black))  ///
legend(label(1 "Debit & Smartphone" ) label(2 "Smartphone")) ///
ytitle("Share of households", height(5)) ///
graphregion(color(white)) bgcolor(white) 
graph export "$dirout/FigureE1c.eps", replace 

use temp.dta, clear
collapse (mean) celular celular_aplicaciones utilizo_debito smart_and_debit [pw=FACTOR], by(trabajo_bins)

gen celular_aplicaciones2=celular_aplicaciones-smart_and_debit
graph bar smart_and_debit celular_aplicaciones2, over(trabajo_bins) stack ///
bar(1, color(cranberry) lcolor(black)) bar(2, color(midblue) lcolor(black))  ///
legend(label(1 "Debit & Smartphone" ) label(2 "Smartphone")) ///
ytitle("Share of households", height(5)) ///
graphregion(color(white)) bgcolor(white) 
graph export "$dirout/FigureE1a.eps", replace 

